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ABSTRACT 

We consider the acceleration of charged particles at the ultra-relativistic shocks, with 
Lorentz factors r s 3> 1 relative to the upstream medium, arising in relativistic fireball 
models of gamma-ray bursts (GRBs). We show that for Fermi- type shock acceleration, 
particles initially isotropic in the upstream medium can gain a factor of order in 
energy in the first shock crossing cycle, but that the energy gain factor for subse- 
quent shock crossing cycles is only of order 2, because for realistic deflection processes 
particles do not have time to re-isotropise upstream before recrossing the shock. 

We evaluate the maximum energy attainable and the efficiency of this process, 
and show that for a GRB fireball expanding into a typical interstellar medium, these 
exclude the production of ultra- high-energy cosmic rays (UHECRs), with energies in 
the range 10 18 - 5 -10 20 5 eV, by the blast wave. We propose, however, that in the context 
of neutron star binaries as the progenitors of GRBs, relativistic ions from the pulsar 
wind bubbles produced by these systems could be accelerated by the blast wave. We 
show that if the known binary pulsars are typical, the maximum energy, efficiency, 
and spectrum in this case can account for the observed population of UHECRs. 
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. ' 1 INTRODUCTION 

X 

5^ ■ The observation of some hundred cosmic ray events to date 
d ' with energy above the Greisen-Zatsepin-Kuz'min cutoff, 
E GZK — 10 19 ' 5 eV, (e.g. Takeda et al. 1998, and references 
therein) has sparked renewed interest in their origin. Ultra- 
high-energy cosmic rays (UHECRs), with energies in the 
range E ~ io 18 ' 5 -10 20 ' 5 eV, are generally believed to be 
extragalactic in origin, based on their harder spectrum, 
isotropic arrival directions on the sky and the fact that they 
are not confined by the Galactic magnetic field. One class of 
models involves continuous production of UHECRs at large 
shocks, such as those associated with powerful radio galax- 
ies, or with ongoing large-scale structure formation in clus- 
ters of galaxies (e.g. Norman, Melrose & Achterberg 1995, 
and references therein). A major difficulty with these mod- 
els is that the number of possible sources within the volume 
contributing to the flux above E GZK , which has a radius 
-Dmax ~ 50-100 Mpc, is too small to explain the observed 
number of independent events. 

An alternative model (Waxman 1995a; Vietri 1995; Mil- 
grom & Usov 1995) considers impulsive production in the 
sources of gamma-ray bursts (GRBs). Observations of X- 
ray and optical afterglows of GRBs (van Paradijs et al. 



1997; Metzger et al. 1997) have confirmed the cosmological 
origin of the phenomenon; the energy associated with each 
event is then thought to be £ GRB ~ 10 51 -10 53 erg. Waxman 
(1995a) and Vietri (1995) noted that assuming comparable 
efficiencies for gamma-ray and UHECR production, the esti- 
mated GRB rate, Qgrb ~ 10 -8 Mpc -3 yr _1 , implies a flux 
of UHECRs reaching Earth from within D max remarkably 
similar to the one observed. Another, perhaps more com- 
pelling, argument is that the dispersion in UHECR arrival 
times due to small-angle deflections in the intergalactic mag- 
netic field implies that at any one time, enough GRB sources 
contribute to the UHECR flux to account for the observed 
number of independent arrival directions (Miralda-Escude 
& Waxman 1996; Achterberg et al., in preparation). 

Relativistic fireball models for GRBs involve an ultra- 
relativistic blast wave with Lorentz factor r s ~ 10 2 -10 3 
bounding the fireball (Rees & Meszaros 1992), and inter- 
nal mildly relativistic shocks (r s ~ 2-10) due to unsteady 
outflows (Rees & Meszaros 1994). The quenching of inter- 
stellar scintillation observed in GRB radio afterglows (Frail 
et al. 1997) confirms the relativistic expansion of the source. 
Vietri (1995) proposed that these relativistic shocks can be 
sites of efficient particle acceleration up to energies exceeding 
10 20 eV. In particular, he argued that an ultra-relativistic 
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shock with Lorentz factor F s will lead to an energy gain per 
crossing cycle Ef/Ei ~ rf, where Ei and Ef are the ini- 
tial and final particle energies. If such an energy gain could 
be obtained repeatedly at the ultra-relativistic blast wave, 
UHECR energies would be reached in only a few cycles. 

In this letter, we first consider the Fermi acceleration 
process at ultra-relativistic shocks in some detail; apart from 
the recent simulations by Bednarz and Ostrowski (1998), 
earlier calculations of such relativistic shock acceleration 
have concentrated on mildly relativistic shocks (r s ^ 10). 
We distinguish between the energy gain in the initial and 
subsequent shock crossing cycles, and show that for phys- 
ically realistic particle deflection processes upstream, the 
latter yield only an energy gain E{/Ei ~ 2. We estimate 
the time scale for the acceleration process, and examine the 
maximum energy attainable, as well as the global energetics, 
for UHECR production scenarios at relativistic blast waves, 
first in a typical interstellar medium, and then in the plau- 
sible environment of a pulsar wind bubble. 



2 ULTRA-RELATIVISTIC SHOCK 
ACCELERATION 

We consider an ultra-relativistic shock of Lorentz factor 
r s ^> 1 relative to the upstream medium. Assuming the 
fluid is weakly magnetised, the shock jump conditions then 
imply that the shock velocity relative to the downstream 
medium reduces to c/3, while the relative Lorentz factor of 
the downstream and upstream media satisfies F r = r s /\/2 
(e.g. Blandford & McKee 1976). 

2.1 Energy gain: initial vs. repeated crossings 

At non-relativistic shocks, the standard scenario of parti- 
cle acceleration (Krymskii 1977; Bell 1978; Axford, Leer 
& Skadron 1978; Blandford & Ostriker 1978) assumes that 
both the upstream and downstream media contain magnetic 
fluctuations which tend to isotropise particles by elastic scat- 
tering in the respective fluid rest frames. 

Assuming that the same general principle operates at 
a relativistic shock, the energy gain per shock crossing cy- 
cle follows from relativistic kinematics (e.g. Peacock 1981; 
Achterberg 1993). For ultra-relativistic particles (of Lorentz 
factor 7 ^> 1), it can be written in terms of initial and final 
energies Ei and Ef as: 

g = r 2 (i-A^„) (i + ft^ u ) = j, (i) 

where the first and second equalities apply to shock cross- 
ing cycles which begin and end in the upstream and down- 
stream media, respectively. Here /3 r is the velocity of the 
downstream medium (in units of the speed of light) rela- 
tive to upstream, r r is the associated Lorentz factor, and 
the quantities /i^d and fi'_, n measure the cosine of the angle 
between the particle velocity and the shock normal, when 
the particle crosses the shock into the downstream and up- 
stream respectively. We use the convention that the shock 
normal points into the upstream medium, so that /Z_, u > 0. 
Primed and unprimed quantities are respectively measured 
in the downstream and upstream rest frames. The only as- 
sumption is that scattering is elastic in the local fluid frame. 



Kinematics require that 1 > //_ u > /3' s = |, so that the 
factor (1 + /3,/i'^u) in (|l|) is of order unity. If fi^d is more 
or less isotropically distributed, as would be the case for 
a population of relativistic particles already present in the 
undisturbed upstream medium, the factor (1— /? r /x_,a) is also 
of order unity, and energy gains Et /Ei ~ T 2 can be achieved, 
as envisioned by Vietri (1995). A similar conclusion holds for 
particles which are initially non-relativistic upstream. 

For all but the initial shock crossing into the down- 
stream medium, however, the distribution in /i^d for a rela- 
tivistic shock will be highly anisotropic (Peacock 1981). For 
an ultra-relativistic particle with Lorentz factor 7 2> r s , 
the Lorentz transformation of the incidence angle reduces to 
fi' « (2-r s 2 2 )/(2 + rf6> 2 ) when 9 = cos" 1 fi < 1. The kine- 
matic condition /iV, u > | is thus equivalent to 9^ u < 1/T B , 
which defines the shock 'loss cone' in the upstream frame. 
We show below that for physically realistic deflection pro- 
cesses upstream, the particle cannot be deflected very far 
beyond this loss cone before the shock overtakes it, so that 
the angle 9^d ~ 1 /r s as well. This severely limits the energy 
gain attainable for all but the initial shock crossing: when 

<§C 1, the downstream energy gain in ([j]) reduces to 

E\ = 1+/V-.U _ 2 + r 2 s et d 
which is of order unity if 9^ d ~ l/r s . 

2.2 Upstream dynamics 

We consider two mechanisms for the upstream deflection of 
particles needed to allow repeated shock crossings: regular 
deflection by a large-scale magnetic field, and scattering by 
small-scale magnetic fluctuations. 

2.2.1 Regular deflection 

We first examine the case where the magnetic field may 
be considered uniform over the region sampled by the par- 
ticle in its excursion upstream. We take the shock veloc- 
ity to define the z-direction, and without loss of general- 
ity assume that the magnetic field lies in the x — z plane, 
B — (B ± , , By), with qB ± > 0. 

The equation of motion for the velocity (3 = v/c of an 
ultra-relativistic particle of charge q can be solved approxi- 
mately by expanding /J z = /i w 1— if? 2 to lowest order in (3 X 
and P v , as 9 2 « /3 2 +/3 2 = 0(1/T 2 ) throughout the particle's 
orbit upstream. Then as long as B± 3> B||/r s , which will 
be the case for all but a fraction ~ 1/T 2 of possible mag- 
netic field orientations, we can ignore the terms in By in the 
equation of motion. If we denote by /3i and /3f the ingress 
and egress particle velocities, we find that (3t x = /3i x and 

Piy = ~ ~ A, + \J-^-Wl-\Piy ■ (3) 
Since (5 2 x + /3? y < 1 /F 2 , eq. (|^) implies 

i<r s f^ d <2 i>^ d >_i. (4) 

The energy gain ratio can then be obtained by substi- 
tuting 9 2 ^ d fa f3 2 x + /3 2 y in ^ . It has its maximum at [3 ix = 
and r s (3 iy w -0.27, where E' { /El w 2.62 . 
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2.2.2 Direction- angle scattering 

We now consider the case where the magnetic field upstream 
is not uniform but irregular. In the simplest such model, the 
field consists of randomly oriented magnetic cells with field 
amplitude B and radius (correlation length) £ c . If the parti- 
cle's gyration radius, r g = E/qB, satisfies r s /T s < £ c , it will 
be deflected sufficiently to recross the shock within a single 
cell, and the regular deflection regime above will apply. In 
the opposite case, r s /T s ^> £ c , the particle's momentum di- 
rection will diffuse in time, with angular diffusion coefficient 
Do = c ^c/(3rg) (Achterberg et al., in preparation). 

In this scattering regime, a particle initially crossing the 
shock at an angle 9^ u will, after a time t, be at an average 
distance upstream relative to the shock given by 



(z(t)} - Zs(t) » ct 



1 

2ff 



Dot 



(5) 



where we have used the fact that 9 = 0(1/T S ) and expanded 
to lowest order in 1/IV Defining the typical upstream resi- 
dence time t u as the solution of (z(t u )) ~ z s (t u ), the typical 
downstream recrossing angle may be estimated as 

2 



9l d ) ~ <*'(*„)) 



f 

' — * U ■ 



(6) 



The typical direction angle upon recrossing the shock, 
<#l d } 1/2 , is thus again of order 1/IV Substituting ^ into 
the energy gain formula (j2|) shows that the energy gain 
typically reaches its largest values when 9^ u = 0, where 
E[/E[ ~ 2. These results are approximate, rather than ex- 
act averages, in that they neglect the correlation between 
individual shock recrossing times and crossing angles. 

2.3 Acceleration time scale 

As the energy gain per shock crossing cycle, AE = Ef — E\, 
is comparable to the initial energy Ei, the acceleration time 
scale will be of order the cycle time, which is the sum of the 
typical upstream and downstream residence times t u and 
t<j. In the case of deflection by a uniform magnetic field, 
the upstream residence time is the time required for the 
particle to be deflected by an angle of order 1/F S [eq. (0)], 
i.e. t u ~ (r^^)- 1 = E/(qV s B ± c). 

The downstream residence time will depend on the 
downstream scattering process, but if we assume Bohm dif- 
fusion, it can be roughly estimated as the gyrotime, i.e. 
t' d ~ (ii^ 1 = E' /(qB'c). If the downstream magnetic field 
is simply due to the compression resulting from the shock 
jump conditions, B' ~ B' ± = v^rsB^, then taking into ac- 
count the Lorentz transformation of the energy for a typical 
particle emerging upstream (E' ~ E/T s ) and the transfor- 
mation of the time interval between two events occurring 
at the shock (t' d ~ td/T s ), one finds that t c \ ~ t u . Moreover, 
turbulence downstream could amplify the magnetic field well 
above the value resulting from simple compression; assum- 
ing it reaches equipartition with the thermal pressure down- 
stream, one can show that in this case B' ~ (c/va)F s B, 
where «a is the upstream Alfven speed. This yields a corre- 
spondingly shorter downstream residence time: 

t d ~ n f „ (7) 



q c^L sV 47r /° 
where p is the upstream mass density. 



The full shock crossing cycle time is thus dominated by 
t u in all cases of interest. In the case of scattering upstream, 
where r g 3> £ c , solving (||) yields a typical upstream resi- 
dence time t u ~ r 2 /(cF 2 £ c ). This may be combined with the 
regular deflection case in the single expression 



tu 



E 



(8) 



qT s Bc 

For a given upstream magnetic field amplitude B, the regu- 
lar deflection regime thus constitutes an absolute lower limit 
on the acceleration time for repeated shock crossings. 



3 FIREBALLS IN THE GENERAL ISM 

We now review the properties of the relativistic blast wave 
driven into the surrounding interstellar medium in fireball 
models of gamma-ray bursts. After an acceleration phase, 
an initially radiation-dominated fireball enters a relativis- 
tic 'free expansion' phase, driving a blast wave with the 
approximately constant Lorentz factor r s ~ v2?7, where 
r\ = £grb /(Mc ), M being the baryonic mass in which the 
fireball energy £grb is initially deposited (Piran, Shemi & 
Narayan 1993; Meszaros, Laguna & Rees 1993). The param- 
eter rj is generally assumed to lie in the range 10 2 -10 3 . 

If cooling of the shocked material can be ignored, this 
is followed by an adiabatic deceleration stage in which 
the blast wave Lorentz factor decreases with radius R s as 
T s oc R~ s/2 (Blandford & McKee 1976). The transition 
between the free expansion and deceleration phases occurs 
roughly at the radius Rd at which the energy in the swept-up 
surrounding material becomes comparable to £grb: 



Rd 



3 £ c 



4n r/ 2 



1/3 



(9) 



where e is the total energy density of the surrounding mate- 
rial; for a typical interstellar medium of number density n, 
this is simply the rest-mass energy e s 
quent evolution of r s is then given by 



_ / r? \ -3/2 / 



3 So 



2tt 



1/2 



R. 



-3/2 



(10) 



and thus once R s > Rd the behaviour of T s is independent 
of the value of r\. 

We now consider various scenarios for the acceleration 
of UHECRs by the relativistic blast wave of a fireball ex- 
panding into the general interstellar medium. 



3.1 Maximum energy for Fermi acceleration 

In order for particles to be accelerated beyond the initial 
boost by shock acceleration of the Fermi type, repeated 
shock crossings must occur. The acceleration time, which 
we argued above is dominated by the upstream residence 
time t u , must therefore be shorter than the age of the blast 
wave measured in the upstream frame, R s /c. Using (g), the 
maximum energy for which this can be the case is 



E <> qBFs x min (^R s , y/Ut, 



(11) 



The most favourable regime energetically is that of regular 
deflection upstream, which requires £ c <; 7? a ; in what follows 
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we will assume that this is the case. Note that the maximum 
energy ( Jj"l| ) is larger by a factor r s than that resulting from 
a simple geometrical comparison of the gyration radius with 
R s . This is due to the fact that a particle typically only 
executes a fraction ~ T" 1 of a Larmor orbit upstream before 
recrossing the shock. 

The behaviour of r s as a function of Rs implies that the 
highest energy in (jll]) is reached at the transition between 
the free expansion and deceleration phases, Rs — Rd- Using 
the definition (^) of R c \, we obtain: 

£^5x 10 15 Z B_ 6 £52 /3 % /3 < 1/3 eV , (12) 

for ions of charge q = Ze, where B-e = B/ (lO _6 G), £52 = 
£grb/ (l0 52 erg), 773 = 77/IO 3 , and n = n/ (lcm~ 3 ). Given 
typical interstellar magnetic fields of a few microgauss and 
the weak dependence on the other parameters, this rules 
out the production of UHECRs by Fermi acceleration at the 
unmodified external blast waves of relativistic fireballs. 



4 FIREBALLS IN PULSAR WIND BUBBLES 

The surrounding medium in which the relativistic fireball 
explodes is probably one modified by the activity of the 
progenitor system prior to the burst event. In the neutron 
star binary merger scenario for GRBs, these progenitors are 
identified with the binary pulsar systems observed in our 
own Galaxy (Narayan, Paczyhski & Piran 1992). The pul- 
sar in these systems can be expected to emit a relativistic 
wind, which over the lifetime of the binary can fill a large 
surrounding volume with shocked relativistic plasma. While 
this plasma will probably be predominantly in the form of 
electron-positron pairs created in the pulsar magnetosphere, 
it has been argued that pulsar winds must also contain ions 
in order to account for the accelerated pair spectrum ob- 
served in the Crab Nebula (Hoshino et al. 1992). The Crab 
Nebula's wisps provide additional observational evidence for 
this ion component (Gallant & Arons 1994). 



3.2 Initial boost of galactic cosmic rays 

3.2.1 Maximum energy 

One remaining possibility for blast wave acceleration of 
UHECRs is that they result from the initial shock crossing 
cycle energy boost of a pre-existing upstream population of 
relativistic particles. This is, in essence, the first of the two 
acceleration mechanisms considered by Vietri (1995). 

This initial boost requires only the time t^ for the parti- 
cle to be scattered downstream, which if equipartition fields 
are generated can be considerably shorter than t u , yielding 
a correspondingly higher maximum energy. Requiring the 
downstream half-cycle time ([?]) for a particle of final energy 
E to be shorter than the age of the fireball, we obtain 

E<,7 x 10 20 Z £5 2 /3 % /3 nJ /6 eV . (13) 

This process can thus attain UHECR energies, provided par- 
ticles with sufficient initial energy to be boosted in this range 
are present upstream. To reach 10 20 eV, this requires rela- 
tivistic particles with energy above ~ 10 14 r/^ 2 eV. 



3.2.2 Global energetics 

The energy invested in boosting pre-existing relativistic par- 
ticles is of order T 2 ecR per unit volume swept up by the blast 
wave, where ecR is the upstream energy density of relativis- 
tic particles. Meanwhile, the blast wave expends an energy 
~ r 2 e per unit volume in shock-heating the surrounding 
medium, which has upstream energy density e w nm p c 2 . 
Thus the fraction of the fireball's energy that can go into 
boosting cosmic rays is of order / ~ e C m/e. Taking the 
interstellar medium values in our Galaxy to be typical, 
ecR ~ leVcm -3 and e ~ 10 9 noeVcm -3 , we see that this 
mechanism can only have a very low efficiency, / ^ 10 -9 . 

The efficient yield of UHECRs required in the GRB hy- 
pothesis, / <; 0.1 (Waxman 1995a; Vietri 1995), could only 
be obtained if a large fraction of the surrounding medium's 
energy density was in relativistic particles. This is in fact 
not implausible in the neutron star binary merger scenario 
for GRBs, as we now argue. 



4.1 Maximum energy for boosted ions 

We will assume the known binary millisecond pulsar sys- 
tems PSR 1913+16, PSR 1534+12, and PSR 2127+11C to 
be typical of GRB progenitors. The spiral-in times of these 
binaries due to gravitational radiation are r ~ 3 x 10 7 yr 
(Narayan et al. 1992), while the periods and period deriva- 
tives of the millisecond pulsars observed in these systems 
yield spindown luminosities ~ 10 33 ergs _1 and characteris- 
tic times ~10 8 yr (Taylor, Manchester & Lyne 1993). Thus 
these pulsars will continuously inject relativistic particles 
into the surrounding medium over the lifetime of the sys- 
tem, with an approximately constant wind luminosity. 

The relativistic flow of the pulsar wind will be ther- 
malised in a termination shock, and this shocked material 
will form a relativistic plasma bubble in the interstellar 
medium. At the ages involved, we expect this pulsar wind 
bubble to be approximately isobaric and in pressure equi- 
librium with the surrounding medium. Once their directed 
kinetic energy is randomised by the shock, the pairs may 
suffer energy losses by synchrotron radiation over time, but 
these are negligible for the ions, which can only suffer adi- 
abatic expansion losses. In an isobaric bubble these losses 
are negligible, so that the ions approximately conserve their 
post-shock energy throughout the bubble. 

The bulk Lorentz factor 7 W of the pulsar wind can be 
estimated by assuming that the entire spindown luminos- 
ity E goes into a wind composed of electron-positron pairs 
and ions, with number fluxes parameterised in terms of the 
Goldreich- Julian flux JV GJ : E = (ftimi + h±rn c )'y v ,c 2 Ngj , 
where mi and m e are the ion and electron masses and hi 
and h± are the ion and pair multiplicities. The wind en- 
ergy per Goldreich-Julian charge is comparable to the polar 
cap potential, and can be estimated as E/N GJ ~ e(_E/c) 1//2 
(see, e.g., Michel 1991). We expect the ions, although a 
minor component of the pulsar wind particles by number 
(hi <C h±), to form the dominant component by mass 
(hiirii <: h±m c ), so that the fraction of the wind energy 
carried by them, £i = [1 + h±m c / (/iim;)] -1 , will be close 
to unity (Hoshino et al. 1992). We will also assume that 
the ions carry the magnetospheric return current, so that 
hi ~ 1/Z (Gallant & Arons 1994). This then determines the 
typical energy, mi7 w c 2 , of the thermalised post-shock ions. 
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One can show that for an initially isotropic distribution 
of relativistic particles upstream, the average energy gain 
(0) for boosted particles satisfies ^§F 2 < (E { /Ei) < §r 2 , 
the exact value depending on the typical angle (/^L u ) of 
those particles recrossing upstream. Taking a factor 2F 2 to 
be representative, the typical energy of the ions boosted in 
the free expansion phase of the blast wave is 



_E max ~ 2rj 2 m i f„c 2 



10 20 % Z & EH 2 eV 



(14) 

where E Xi = E/ (lO 33 ergs -1 ). UHECR energies are thus 
naturally obtained is this scenario, for typical values of the 
pulsar wind parameters, as long as 77 i> 10 3 . 

4.2 Spectrum and efficiency 

Although the upstream ions have a relatively narrow en- 
ergy spectrum centered around rriij v c 2 , the boosted ions 
will have a broader distribution in energies provided the 
blast wave decelerates within the pulsar wind bubble, so 
that the boosting factor F 2 is not constant. The decelera- 
tion radius Rd is given by (^), with e the energy density of 
the relativistic plasma, while the pulsar wind bubble radius 
Rb is obtained by equating the volume integral of e with the 
total energy output of the pulsar, ft follows that 



77- - 0-3£ B yii , 
tib 



1/3^-2/3 p-l/3 -1/3 
^33 '7 ' 



(15) 



where r = T7 x 10 7 yr is the age of the system. Thus for typ- 
ical parameters, the relativistic fireball will decelerate well 
within the pulsar wind bubble. 

Since E oc F 2 oc RT 3 in the deceleration phase, and 
the number of ions swept up by the blast wave scales as 
diV oc R 2 dR B , the spectrum of the boosted ions will be 

This is consistent with the observed UHECR spectrum 
(Waxman 1995b). It can be shown that the same spectral in- 
dex follows more generally if the ion density is not uniform, 
as assumed here, but obeys a power law in radius. 

The lower bound of the spectrum ( p"ri| ) follows from the 
boosting factor F 2 of the blast wave when it reaches the edge 
of the pulsar wind bubble at Rb, and is 



-En 



3 x 10 18 ZiiEsaE^' rf 1 eV 



(17) 



Thus not only can a large fraction of the fireball energy go 
into boosted ions, but these ions have a power-law spectrum 
extending more or less exactly over the energy range where 
the extragalactic UHECR component is observed. 



5 SUMMARY 

In this letter, we examined Fermi-type acceleration at an 
ultra-relativistic shock, and found that particles with ini- 
tially isotropic momenta upstream can increase their energy 
by a factor of order F 2 in the initial shock crossing cycle. 
In all subsequent shock crossing cycles, however, we showed 
that the particle energy typically only doubles. This is due 
to the fact that particles do not have time to re-isotropise 
upstream before being overtaken by the shock, which we 
demonstrated in the specific cases of a large-scale ordered 
magnetic field and of small-scale magnetic fluctuations. 



We argued that the maximum energy that can be 
reached by repeated shock crossings at the unmodified rel- 
ativistic blast wave from a GRB fireball is well below the 
UHECR range. Pre-existing relativistic particles of sufficient 
energy can nonetheless be boosted to UHECR energies in 
the first shock crossing cycle. For a fireball expanding into a 
typical interstellar medium, where galactic cosmic rays pro- 
vide the seed particles, we showed, however, that this process 
is too inefficient to account for UHECR production. 

We proposed that the blast wave instead expands into 
a pulsar wind bubble produced by the progenitor system, 
a plausible hypothesis in the neutron star binary merger 
scenario for GRBs. We showed that for parameters typical 
of the neutron star binary systems observed in our Galaxy, 
relativistic ions in the pulsar wind bubble can be efficiently 
boosted by the blast wave to energies exceeding 10 20 eV. We 
argued that these boosted ions would have an E~ 2 spectrum, 
extending down in energy to about lO 18 5 eV. 
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